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NASA is concerned with protecting astronauts from the effects of galactic cosmic radiation and has expended substantial effort in 
the development of computer models to predict the shielding obtained from various materials. However, these models were only 
developed for shields up to about 120 g/cm 2 in thickness and have predicted that shields of this thickness are insufficient to provide 
adequate protection for extended deep space flights. Consequently, effort is underway to extend the range of these models to thicker 
shields and experimental data is required to help confirm the resulting code. In this paper empirically obtained effective dose 
measurements from aircraft flights in the atmosphere are used to obtain the radiation shielding function of the earth’s atmosphere, a 
very thick shield. Obtaining this result required solving an inverse problem and the method for solving it is presented. The results are 
shown to be in agreement with current code in the ranges where they overlap. These results are then checked and used to predict the 
radiation dosage under thick shields such as planetary regolith and the atmosphere of Venus. 
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I. Introduction 

Even before Apollo, the dangers of deep-space radiation were recognized, and significant effort has since been spent in 
quantifying this radiation, in understanding its effects on humans, and in finding ways to ameliorate it — and we now know that 
the situation is more perilous than previously thought. There are still unknowns, but the recent report (2012) from the National 
Academies on NASA Space Technology Roadmaps and Priorities summarizes the situation: “...models predict that crewed 
missions beyond low Earth orbit (LEO) would be limited to three months or less because of adverse health impacts.... ”[1] Yet 
three months is insufficient for manned space flight beyond LEO, except for lunar sorties, causing deep-space missions to be 
presently unachievable. Radiation effects pose such a fundamental obstacle to man’s exploration of the solar system that the 
National Academies state that “Space Radiation Health Effects” (including investigating protection technologies) should be 
NASA’s second highest technical challenge (behind “Improving Access to Space”). 

The majority of this health detriment is caused by the accumulated effect of galactic cosmic radiation (GCR). GCR is composed 
of high-energy nuclei ranging from protons to uranium that enter the solar system from all directions. When these “projectiles” 
encounter matter, either in the structure of a spacecraft or within a human body, they collide with and shatter nuclei, causing a 
shower of secondaries composed of charged particles, neutrons, and gamma rays — all of which are damaging to human tissue. 
Neurons are harmed irreparably, cataracts can start to develop, and the risk of dying from cancer, i.e. the risk of exposure- 
induced death (REID) rises [2]. 

With this picture in mind, researchers in the 1990s began to extol the benefits of hydrogen-containing materials as shields against 
GCR, because hydrogen has no neutrons to add to the secondaries and it is the most effective material per unit mass at slowing 
down ions [3]. For a few years it appeared that these materials might provide adequate shielding, and projects were funded to 
study high-hydrogen-containing materials, especially polyethylene [4]. But the models of the early 2000s did not adequately 
account for the interaction of radiation with the human body and its tissues. Refinements in these models led to predictions of 
high levels of radiation [2,5,6] and are the basis for the statement quoted by the National Academy in the introductory paragraph 
above. It was shown that adding substantial polyethylene could actually increase the effective radiation dose received by an 
astronaut and that 120 g/cm 2 of polyethylene (more than a meter of material) could not provide the shielding necessary to allow a 
mission to Mars [see Fig. 6.5 in ref 2]. 

In the face of this picture a reasonable question to ask would be, “How much shielding is sufficient to protect an astronaut for a 
given mission?” However, even though the various radiation codes are highly developed [7] they have only been verified against 
in-flight and accelerator data out to about 120 gm/cm 2 and beyond this the models use exponential extrapolations. Model 
verification data to high shielding mass is needed to have confidence for spacecraft and surface habitat design. One of us 
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(Koontz) recognized that such data might be obtained by manipulating information from the Federal Aviation Authority (FAA). 
The FAA has developed a CARI-6 computer program that supplies GCR effective dose values for any location in the world from 
ground level up to 60000 feet [8,9] and is available as a free download. The data within the CARI-6 is a combination of transport 
code, human body modeling, and dosimeter measurements [10] and represents real world measurements of the effective dose an 
individual receives from galactic cosmic radiation when protected by a thick atmospheric shield. The problem is then to extract 
shielding functions for the earth’s atmosphere from the CARI-6 data, requiring solving an interesting inverse problem. 

II. Theoretical Analysis 

The data supplied by the CARI-6 program describes the effective radiation dose received by a human from galactic cosmic 
radiation, which is isotropic, arriving at the earth equally from all directions. Yet, to an observer at some height h above ground 
level, this total effective radiation, /?(//), is the sum of radiation that is not spatially uniform, having passed through different 
mass lengths of the earth’s atmosphere. Mathematically this can be expressed as 


In nil 

R(h) = | J S(m(0,<p,h))sin(0)d0d<p (1) 

0 0 

where S(m) is the effective radiation dosage per steradian received by the observer after having passed through an atmospheric 
mass, m(0,</>,h ) , (units of mass/area) where 6 is the elevation angle, measured from the zenith as shown in Figure 1, and <p is 
the azimuthal angle. In the discussion below 4/rS(m) will be referred to as the atmospheric radiation shielding function since this 
corresponds to the total effective radiation dose seen when surrounded by a spherically symmetric shell of atmosphere. 

The elevation angle in Equation 1 is only integrated to n / 2 radians, corresponding to the horizon and accounting for the 
radiation shielding provided by the planet itself. Also, in all cases of interest in this paper the atmospheric mass will be assumed 
to be independent of azimuth, allowing equation (1) to be rewritten as 


/r/2 

R(h) = In J S(m(0, h))sm(0)dG . (2) 

0 

Mathematically, the problem here is to solve for the function S(m) assuming R(h) is known. Problems such as this are referred 
to as inverse problems and commonly arise in 3-D imaging problems, such as tomography [11] and sonar [12], but can appear in 
diverse fields. A common approach to solving these problems, which we will develop below, is to find a way to rewrite the 
integral as a matrix multiplication and then calculate the inverse matrix. 

The atmospheric mass function, m(0,h) , represents the amount of atmospheric mass between an observer located at a height h 
above the ground, and a point above the atmosphere when looking at an angle 0 down from the zenith; along the s direction 
shown in Figure 1. In order to obtain an expression for this function we will assume an isothermal atmosphere in which pressure 
and density decay exponentially with increasing height. The analytic expressions in this model are given by 

P(h) = P(0)exp[-h /H] and p(h) = -P\h) / g = (P(0) / ( gH))exp[-h / H) = p(0)exp[-h / H ] (3) 

where P(h) is the pressure as a function of height, p(h) is the atmospheric density as a function of height, g is the 
gravitational acceleration (9.8 m/sec 2 ), and H is the scale height. The scale height of the atmosphere is given by H = kT / (Mg) 

where k is Boltzmann’s constant ( 1.38x1 0" 23 J/K), T is the average atmospheric temperature in Kelvin (we use 260 K as an 
average atmospheric temperature), and M is the mean molecular mass of dry air (equal to the effective number of atomic mass 
units in the atmosphere (29) times 1.66x1 O' 27 kg). These values yield a scale height of 7600 meters. In addition we will assume 
that the pressure at ground level, P(0) is given by 101 ,000 Pascals. 



3 



Figure 1. This figure defines the variables used in the calculation of the 
atmospheric mass function. 

Now, referring to Figure 1 , the mass function is given by the integral of the atmospheric density along the line of sight direction 
shown by the s parameter. Writing this explicitly and incorporating the density function in Equation (3) yields 


m(6, h) = J p(s)ds =J p( 0) exp [-z(0,s) / H] ds . (4) 

0 0 

However, the density is a function of height, z , and the integral is over the line of sight variable, s , so we need to determine the 
function z(0,s ) . This can be done by using the law of cosines and noting that all upward directed vectors pass through the center 
of the earth, 

z(s,0) = +Jh 2 +2hr e +/; 2 +s 2 + 2(h + /;)scos(0) . (5) 

where r, is the radius of the earth, 6.378xl0 6 meters. Combining equations (4) and (5) yields the final expression for the 
atmospheric mass function 


m(0,h) = p(0)Jexp[|/; -yjh 2 +2hr e +/; 2 +s 2 +2(/? + /;)scos(0) j / H]ds (6) 

o 

Figure 2 shows three plots of the mass function, for heights of 0, 9000 meters, and 1 8000 meters above the ground versus angle 
from vertical. As the angle approaches the horizon at n 1 2 radians the mass increases rapidly as expected. Note that this plot is 
in units of g/cnr, rather than the SI kg/nr, because most of the radiation shielding literature uses CGS rather than MKS units. 



4 



Figure 2. Plots of the total atmospheric mass between an observer and space versus the angle in 
radians measured down from vertical. The highest plot corresponds to ground level, the middle plot 
is at a height of 9000 m, and the lowest plot is at a height of 18000 m. 


Having derived an expression for the atmospheric mass function, the inverse problem given in Equation (2) can now be 
converted into a matrix expression. Start by selecting a set of discrete masses, m l , based on the masses provided by Equation (6), 
such that the radiation shielding function, S(m ) , can be approximated as being constant over each mass interval. Let m, be the 
smallest mass and m n+] be the largest so that there are n intervals and the radiation shielding function, S(m) can be 
approximated as 5(m y ) for each interval between masses m t and w (+ , . Now chose a discrete height, h j , and for each discrete 
atmospheric mass, m i , at that height let G t be the angle that corresponds to that mass. So between each angle interval, 0. to 0 i+l , 
the atmospheric mass, m . , seen by an observer looking at that angle is approximately constant. In addition, if at some height the 
masses are too small to occur then set the corresponding 6 t = 0 and if the masses are too large then set the 
corresponding 0 % = n / 2 . Doing this allows Equation (2) to be rewritten as 

R(h J ) = "sin(O)d0 = 2n'£S(m l )(cos(0 jJ+l )-cos(0 j j)) = Y t S(m i )Q J j (7) 

1 = 1 ' i=\ 1=1 

where the matrix 0 y . describes the angular ranges at each height, , through which an observer sees radiation corresponding to 
S(m t ) . Writing this in matrix notation yields 


R = SS (8) 

which can be solved to obtain the atmospheric radiation shielding function, AnS (m ) , but care is required. Simply inverting the 
matrix, Qj . , in Equation (8) can yield chaotic results because small noise or errors in the input data, R(hj), can produce large 
variations in the output data, S(m t ) . The present case is no exception and its explicit solution will be discussed in the next 
section. 
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III. CARI-6 Data Analysis 


The FAA’s CARI-6 computer program supplies the effective dose values for any location in the world, from ground level up to 
60000 feet, averaged over any month going back to 1958. For the purposes of the present study we chose our location to be the 
north geomagnetic pole (80 degrees N, 72 degrees W) in order to minimize the impact of the earth’s magnetic field on the 
propagation of the GCR. We then chose two months, June of 2000, corresponding to a solar maximum where the sun’s magnetic 
field was strong, reducing the amount of GCR that reaches the earth, and June of 2009, corresponding to a solar minimum where 
the GCR impact on the earth was high. We tabulated this data for every 2500 feet from 2500 feet to 50000 feet (we dropped the 
ground point because it appeared to have a background radiation contribution) and then every 1000 feet from 50000 feet up to 
60000 feet. The interval was decreased at the higher elevations to help reduce uncertainty in the lower mass range. This data is 
shown in Figure 3 with the lower plot corresponding to solar maximum and the upper plot corresponding to solar minimum. 




0 5000 10000 15000 


height(m) 

Figure 3. This is the CARI-6 effective radiation doses from GCR versus height at the 
north geomagnetic pole for June 2000, a solar maximum when the radiation levels 
were low, and for June 2009, a solar minimum when the radiation levels were high. 


A polynomial fit to the CARI-6 data was found in order to supply effective radiation doses at arbitrary heights. Then, using the 
results of the theory section, a discrete matrix form for the inverse problem was developed allowing the inverse problem to be 
solved using the Mathematica routine “LeastSquares”. This specifies that the square root of the sum of the squares of the 
differences between the fit and the data be used as a norm to find the best solution to the inverse matrix problem. Doing this 
yielded a pair of atmospheric radiation shielding functions, one for solar maximum and one for solar minimum. These two 
functions are the primary results of this effort and are shown in Figures 4 and 5. 

Figures 4 and 5 show the effective radiation dose that an astronaut would receive in deep space when surrounded by a shell of 
atmosphere with a given mass density. By deep space we mean a location well outside of the earth’s protective atmosphere and 
magnetic field, yet at a distance from the sun similar to that of the Earth. This form is typical of that found in the literature and 
allows a direct comparison with previously published results. For example Cucinotta [13-Fig. 6.6] shows the effective dose for a 
male behind an aluminum shield with masses ranging from 0 to 120 g/cnr. For solar maximum his results show an effective 
radiation dose between 270 mSv/year and 320 mSv/year (depending on the model used) at 120 g/cnr shielding. Our solar 
maximum plot is consistent with this result showing an effective radiation dosage of 276 mSv/year at 120 gm/cm 2 . However, 
our results extend beyond 120 g/cm 2 showing a peak radiation dose in the 160-180 g/cm 2 range of 287 mSv/year. For solar 




minimum Cucinotta shows a radiation dosage of between 430 mSv/year and 520 mSv/year at 120 g/cm 2 mass. Our results are 
larger, see Figure 5, showing about 580 mSv/year at 120 g/cm 2 . This reason for this discrepancy is unclear, but might be due to 
differing choices of solar minimum. 



Mass (g/cm A 2) 

Figure 4. This plot shows the atmospheric radiation shielding function at solar maximum, i.e. this 
is the radiation dose an astronaut would receive in deep space at solar maximum when surrounded 
by a shell of atmosphere of given mass thickness. 



Mass (g/cm A 2) 

Figure 5. This plot shows the atmospheric radiation shielding function at solar minimum, i.e. 
this is the radiation dose an astronaut would receive in deep space at solar minimum when 
surrounded by a shell of atmosphere of given mass thickness. 
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The unique feature of the atmospheric radiation shielding functions shown in Figures 4 and 5 is that they provide data on the 
high mass shielding of the atmosphere. For atmospheric masses larger than 200 g/cnr the curves are roughly exponential, 
decaying about a factor of 3 for every 150 g/cm 2 of additional mass. Also, these results are applicable to shielding materials other 
than the constituents of the earth’s atmosphere. Studies have shown that most non-hydrogen containing materials have similar 
shielding properties [see. Figure 4 in reference 5] when compared by mass density. Consequently, the atmospheric radiation 
shielding functions in Figures 4 and 5 can be used with a reasonable degree of confidence to predict the shielding effectiveness 
of other materials. For example, a question commonly raised by the space community is, “How much regolith is needed on the 
Moon to protect astronauts for long time periods?” These plots show that 1000 g/cm 2 is more than sufficient to provide long 
term radiation protection. So, if one intended to use the lunar regolith as a radiation shield and its density is 3 g/cnr then a layer 
between 3 and 4 meters thick would provide maximal protection with no need to dig deeper. 

NASA is currently considering setting 150 mSv as a maximum allowable lifetime dose for an astronaut and these curves indicate 
the amount of shielding that would be necessary to achieve that for a given time period in deep space. The resulting amounts of 
material may be excessive for use as shielding on a spacecraft, though a recent publication has considered this possibility [14]. 


IV. Checking the Result and Predicting the Radiation Dosage in the Venusian atmosphere 

Before considering other planetary atmospheres, we begin with Earth in order to check our results and to test the forward 
calculation process. Inserting the atmospheric radiation shielding functions into Equation 2, the radiation dosage seen at various 
heights over the north geomagnetic pole can be calculated. The results are shown in Figure 6 along with the CARI-6 data, 
demonstrating that the shielding plots found through the inverse method described above do indeed yield the correct effective 
radiation dosages when used in the forward problem. 



0 5000 10000 15000 

Height (m) 


Figure 6. This plot shows the effective radiation dosage for heights above the geomagnetic north pole on 
Earth, for solar minimum (the upper plot) and for solar maximum (the lower plot). The points are the 
CARI-6 data. This plot shows that the solutions to the inverse problem can be used in the forward 
direction to recover the original data. 
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Now that the forward process has been verified for Earth’s atmosphere we can proceed to model the total effective radiation 
dosage that an astronaut would receive at various heights in the atmosphere of Venus. We chose Venus to study because other 
entities within the solar system either have no appreciable atmosphere (e.g. Mercury, Mars) or have primarily hydrogen 
atmospheres (e.g. Jupiter) for which the above analysis is not applicable or have an atmosphere but are at a much larger distance 
from the sun where the GCR radiation levels are not known (e.g. Titan). For clarification, Mars has a thin atmosphere that might 
provide some shielding. We calculated the effective radiation dosage for an explorer on the surface using Cucinotta’s results for 
angles where the atmospheric mass is less than 100 gm/cm 2 and our results for angles (near the horizon) where the total 
atmospheric mass is thick. Our analysis showed that the radiation dosage on the surface of Mars is roughly half of that in deep 
space, due to the planet itself blocking half of the sky and the atmosphere providing essentially no shielding (we did not 
calculate the additional radiation caused by GCR interaction with the ground). So for the purposes of radiation shielding the 
Martian atmosphere is ineffective. 

Table 1 shows the pertinent data needed to perform atmospheric effective radiation dosage calculations for Venus where Earth 
and Mars are included for comparison. The average atmospheric temperature, the mean molecular mass, and the gravitational 
acceleration are needed in order to calculate the scale height of the respective atmosphere. Using this with the surface pressure 
allows the surface density to be calculated and then using the body radius in Equation (6) allows the shielding mass to be found. 
Equation 2 can then be used to determine the total effective radiation dosage versus height within that atmosphere. 



Earth 

Mars 

Venus 

Average Atmos. Temp K 

260 

210 

300 * 

Mean Molecular Mass (AMU) 

29 

43.3 

43.4 

Gravitational acceleration (m/s 2 ) 

9.81 

3.7 

8.87 

Scale height (m) 

7600 

10900 

6500* 

Radius (m) 

6378000 

3396000 

6051000 

Surface pressure (Pa) 

101000 

636 

107700 * 

Surface density (kg/m A 3) 

1.35 

0.016 

1.87* 

Notes/issues - 


Lowest point is -7 
km. 

* 50 km above the 
actual surface 


Table 1. The parameters required to derive the Galactic Cosmic Radiation shielding characteristics of 
various atmospheres. 


Venus is sometimes called Earth’s sister planet and this relationship applies to the two atmospheres, assuming the comparison 
with Earth is done 50-55 kilometers above the surface of Venus. The surface of Venus is over 700 K and is extremely hostile, 
but at a height of about 50-55 kilometers the Venusian atmosphere has similar pressure and temperature to the Earth’s 
atmosphere at ground level. The benign condition at this location has caused at least one researcher [15] to suggest colonizing 
Venus, referring to it as the most Earthlike environment in the solar system. So it is not surprising that the effective radiation 
dose starting 50 km above the surface of Venus point is very similar to Earth’s — compare Venus’s effective radiation dose curve 
shown in Figure 8 to Earth’s shown in Figure 6. A floating manned science platform on Venus could range between 50 and 57 
km above the surface for years without exceeding the lifetime radiation dose suggested by NASA. 
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Height (m) 

Figure 8. These are plots of the effective radiation dose within the atmosphere of Venus versus height 
(starting at 50 kilometers above the surface) for solar maximum (lower plot) and solar minimum (upper plot). 


V. Conclusion 

We have shown that by using effective radiation dosage data from the FAA’s CARI-6 program that the atmospheric shielding of 
the Earth’s atmosphere can found. To our knowledge this is the first time that an effective radiation shielding curve has been 
found for very thick, up to 1000 g/cnr, mass densities. This information can provide verification for computer models used to 
predict material shielding properties. In addition, it supplies useful insight into how much material on the moon or Mars would 
be required to construct an effective radiation shield and how much shielding would be required to protect an astronaut in deep 
space. It also allows modeling of the effective radiation dosages expected within the atmosphere of Mars and Venus. As might 
be expected, the Martian atmosphere provides very little radiation protection; however, there is a benign region in Venus’ 
atmosphere, ranging between 50 km and 57 km above the ground where there is a high degree of radiation protection. 
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